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We study the morphology of magnetic domain growth 
in disordered three dimensional magnets. The disordered 
magnetic material is described within the random-field Ising 
model with a Gaussian distribution of local fields with width 
A. Growth is driven by a uniform applied magnetic field, 
whose value is kept equal to the critical value H C (A) for 
the onset of steady motion. Two growth regimes are clearly 
identified. For low A the growing domain is compact, with 
a self-affine external interface. For large A a self-similar 
percolation-like morphology is obtained. A multi-critical 
point at (A c , H C (A C )) separates the two types of growth. We 
extract the critical exponents near A c using finite-size scaling 
of different morphological attributes of the external domain 
interface. We conjecture that the critical disorder width also 
corresponds to a maximum in H C (A). 

I. INTRODUCTION 

Many physical processes involve motion of the interface 
between two domains through a disordered media. Clas- 
sic examples are magnetic domain growth, fluid invasion 
of porous media and fluid segregation in gels .Em Previous 
studies of such systems have shown that changing the 
strength of disorder leads to critical transitions between 
different growth morphologies Hllil Most theoretical work 
on the subject has been done on models with an under- 
lying crystalline lattice of spins or pores. The crystalline 
anisotropy may lead to a faceted growth regime in the 
low disorder limit. As disorder increases, there may be 
a transition to compact growth with a self-affine inter- 
face. In the high disorder limit, growing domains are 
self-similar fractals, characteristic of percolation. 

Studies of driven interfaces in the zero-temperature 
random-field (and random-bond) Ising model (RFIM) 
indicate that the dimensionality, coordination number, 
and distribution of random fields are all important in de- 
termining the sequence of morphological transitions .ETlI 
In two dimensions (2D), there is a transition from self- 
similar to faceted growth at a critical value of disorder, 
if the the distribution of random fields is bounded!] The 
critical behavior is not universal, and can be related to 
the analytic form of the tails in the distribution!] For 
an unbounded Gaussian distribution, percolative growth 
occurs for any finite strength of disorder. In the three di- 
mensional (3D) RFIM with a bounded distribution of dis- 



order all three types of growth morphology, faceted, self- 
affine, and percolative, are found as the strength of disor- 
der increases.^ Analysis of growth probabilities suggests 
that the transition from percolative to self-affine growth 
might be universal, but that the self-affine to faceted 
transition is non-universal and could be suppressed by 
an unbounded distribution of fields. 

The suggestion that faceted growth may be eliminated 
by unbounded distributions of disorder in 3D is sup- 
ported by renormalizatkm group studies of equilibrium 
interface conformations]!!] Calculations for a Gaussian 
distribution of random fields show that faceted interfaces 
are suppressed by any amount of disorder in dimensions 
d < 3. Porous media and some magnetic systems do not 
have any underlying crystalline structure, and it is only 
introduced in the models for computational convenience. 
Thus suppression of faceted growth by unbounded dis- 
tributions of disorder leads lattice models to provide a 
more accurate description of real systems. 

In this paper we explore the effect of unbounded Gaus- 
sian disorder on growth morphology transitions in the 
3D RFIM at zero temperature. We find that this un- 
bounded distribution does indeed eliminate the faceted 
growth regime. There is a single transition from self- 
affine growth at low disorder to percolative growth at 
high disorder. For each value of the disorder we find the 
critical field H c needed to initiate steady growth. The 
critical behavior at the onset of growth is analyzed, and 
the critical exponents are consistent with previous results 
for bounded distributions of random fields. 

The transition between self-similar and self-affine 
growth is shown to be a multi-critical point. We identify 
lengths that diverge as the strength of disorder is varied 
along the line of critical fields, and evaluate critical expo- 
nents. Our analysis reveals problems with previous work 
on bounded distributions of random fields. These studies 
identified a "fingerwidth" with the correlation length that 
diverges at the multi-critical pointE However, our work 
with larger systems shows that this fingerwidth does not 
diverge, and its saturation leads to errors in the deter- 
mination of critical properties. New work will be needed 
to determine whether the self-affine to percolative transi- 
tions for bounded and unbounded distributions are in the 
same universality class. However, we can compare our re- 
sults to studies of magnetic hysteresis in the 311 RFIM 
with Gaussian random fields by Perkovic et aLljii!] These 
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studies use a different growth algorithm, and we find that 
this leads to changes in the critical disorder and critical 
exponents. 

The paper is organized as follows. In Sec. |lj we de- 
scribe the growth model. Studies of the critical field and 
critical expon ents at the onset of motion are presented 
in Sec. 
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T he tra nsition in growth morphology is 
examined in Sec. [II B, and Sec. IV contains our general 
summary and conclusions. 



II. GROWTH MODEL AND ALGORITHM 

The energy of a system of Ising spins (sj = ±1) at the 
sites i of a simple cubic lattice is written as 



H 



SiSj - )Xrji + H)si 



(1) 



<hj> 



The first term on the right-hand side of Eq. (0) repre- 
sents the ferromagnetic coupling between nearest neigh- 
bor spins, with the exchange coupling taken as the energy 
unit. The second term gives the interaction of each spin 
with the uniform external magnetic field H and with the 
random local field rji. The local fields are uncorrelated, 
following a Gaussian distribution function whose width 
A quantifies the degree of disorder: 



P{rj) = (2ttA 2 



-i/2 e - n 2 /(2A 2 ) 



(2) 



As in previous studies,! the simulation cell is a cube 
of side L, with the lattice constant taken as the unit of 
length. Periodic boundary conditions are imposed along 
the x and y directions. All spins are initially anti-parallel 
to the external field, i.e. Si = — 1, except those in the 
bottom layer, z = 1. This layer constitutes the "seed" 
for growth of the +1 domain, which is driven by the 
external field. The orientation of this seed plane does not 
affect critical behavior in the sglikaffine and self-similar 
growth regimes of interest here,QE3tj but does affect the 
faceted arpaith seen for weak, bounded distributions of 
disorder .QiIjIIj 

Growth proceeds through single spin flips at zero tem- 
perature. Whenever a spin flip from Sj = — 1 to Sj = +1 
lowers the total energy, it is implemented. However, only 
spins on the growing interface are allowed to flip. Tl 
restriction is motivated by the fluid invasion problem,! 
and differentiates our .dynamics from the model consid- 
ered by Perkovic et alwEi in studies of hysteresis in mag- 
netic systems. 

We have developed a new memory-efficient algorithm 
that allows us to consider cells with L as large as 1152 in 
the present study. The simulation cell is subdivided into 
smaller cubic cells of side s. The s 3 spins of a sub-cell re- 
main active in memory only while at least one of them is a 
"flippable" interface spin. Growth occurs at fixed driving 
field H . This allows the local random field to be encoded 



into a single byte that gives the minimum number of addi- 
tional neighbors needed to flip the spin at that site. This 
number is decreased by one every time the spin acquires 
a new +1 neighbor. When enough neighbors are pr esent , 
the spin is flipped. For reasons discussed in Sec. IIIB , 
we flipped all spins that were completely surrounded by 
+1 neighbors, regardless of their local field. Such spins 
can not affect growth at any other site, and being able to 
remove the cells containing them from memory allowed 
us to treat larger systems. Preliminary runs where these 
spins were not flipped gave equivalent results. 

At each L and A, domains were grown in an ensem- 
ble of samples with different configurations of {r/i}. In 
each sample, growth was stopped when all spins on the 
interface were stable, or when the interface first reached 
the top of the system (z — L). The morphology of the 
domain was then analyzed as described below. 



III. RESULTS 
A. The critical field 

The external field H provides a driving force that 
causes the domain of +1 spins to grow. If H is too 
small, the domain wall will remain pinned near the bot- 
tom of the system, while large H will cause all spins to 
flip. Figure [l] shows the probability, Ptop, for a domain 
to grow to the top of the cell as a function of H at sev- 
eral different system sizes. Results for two values of A 
corresponding to self-aflme (A = 2.1) and self-similar 
(A = 3.6) growth are shown. For each A, the range of 
H over which the probability rises from to 1 becomes 
narrower as L increases. All of the curves intersect at a 
critical field H C (A) and probability P C (A). If H > H c , 
the probability approaches unity in the thermodynamic 
limit [L — » oo), while the probability vanishes in this 
limit if H < H c . 

PrpvkHigrStudies with a uniform distribution of random 
fieldsoQO'tS have shown that there is a diverging corre- 
lation length £jy ~ \H — H C \- VH as \H - H c \ -0). The 
critical exponent in the self-similar growthxegimeS is con- 
sistent with the 3D percolation exponent Jl3 vjj = ± 
0.02, while v H = 0.75 ± 0.02 for self-affine growth.tH 

Near H c , the probability for the domain to span the 
system should only depend on the ratio of system size 
to £h- This suggests that when probability is plotted 
against (H — H c )L 1 ' l ' H the results for all system sizes 
should collapse onto a universal curve. Fig. || verifies 
that the data from Fig. |l| are consistent with this ansatz, 
and with the values of vh that were found for uniform 
distributions of random fields. 
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FIG. 1. Probability P top that the domain grows to the top 
of the system as a function of H at the indicated system sizes 
L for (a) A = 2.1 (self-afhne) and (b) A = 3.6 (self-similar). 
Four hundred different systems were used for each L and A, so 
the statistical uncertainty is less than ±0.025 (symbol size). 



1 

0.8 
0.6 
0.4 
0.2 

1 



0.2 



_l 1 1 I 1 1 I 1 1 1 


—sax x h— 

x H : 

X 




n — 

x i 


L a 




: e 




: *■ 

_ x& 


A = 2.1 J 


: 

-atx a«xB " 


(a) j 
1 1 1 1 1 f 


J 1 1 1 1 1 1 1 


** J 


- 4 

X 

Z □ 

r ^ 

X 

x>Ox 




A = 3.6 J 
(b) j 



-10 -5 5 

(H-HJL 1 /"- 



10 



FIG. 2. Probability that the domain grows to the top of 
the system as a function of (H — H C )L}I VH for (a) A = 2.1 
and (b) A = 3.6. In (a) H c = 1.4905 and v H = 0.75. For (b) 
H c = 1.335 and vh = 0.88. The symbols for each L are the 
same as in Fig. [j]. Four hundred different systems were used 
for each L and A, so the statistical uncertainty is less than 
±0.025 (symbol size). 



Examination of the spanning probability for other val- 
ues of A shows that -P C (A) is always near 2/3 for self- 
similar ^growth and decreases with A in the self-afhne 
regime.EjTo determine H C (A) we worked with the largest 
accessible system size (L =768 or 1152) and found the 
value of H that gave a spanning probability between 0.5 
and 0.7. This determines H c with an accuracy of about 
0.0001. We found that even larger errors in the growth 
held did not change the morphology of spanning domains 
that is analyzed in the next section. 

In Figure ^| we show the variation of H c with disorder 
for the range of interest in the present study. At small A 
the value of H c increases monotonically. There is a max- 
imum near A = 2.5, and then H c drops monotonically, 
becoming negative for A > 7. 
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FIG. 3. Variation of the critical field H c with Gaussian dis- 
order A. A filled circle indicates the value of A c = 2.52 ±0.03 
determined below. Self-afnne growth occurs for A < A c and 
self-similar growth occurs at larger A. 

B. Transition in Growth Morphology 

As expected, the faceted growth jjegime seen for 
bounded distributions of random heldscl was suppressed 
by Gaussian disorder. The only transition that we ob- 
served was from self-similar growth at high disorder to 
self-afhne growth at low disorder. 

Self-similar growth is isotropic, while self-afhne growth 
has a well-dehned direction. a± long-length scales. Previ- 
ous work on other modelstjlj shows that the transition 
between these two growth regimes is a multi-critical point 
at some H c and A c . An order parameter can be defined 
in analogy to equilibrium magnetic transitions asj-fihe av- 
erage of the unit vector normal to the interface.Qil3 For 
A < A c this average is finite, while for A > A c the 
order parameter vanishes. A correlation length, £, that 
diverges at A c can also be defined. As A approaches A c 
from above (self-similar regime), longer and longer seg- 
ments of the domain wall advance in the same direction. 
In the self-affine regime, deviations from the mean di- 
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rection occur over longer and longer length scales as A 
increases to A c . In the following subsections we exam- 
ine the morphology of domains and use finite-size scaling 
to determine A c and the exponent v that describes the 
divergence of £ as A — > A c . 



1. Fingerwidths 

Previous experimental and theoreticalETEl studies 
have used a simple measure of the range of correlations. 
A fingerwidth w was calculated by examining lines of ad- 
jacent nearest-neighbor spins, and averaging the length 
of contiguous segments of +1 spins (or fluid- invaded 
rcgionsHatrLJ). Results obtained from our simulations 
are given in Figure ||(a). As in previous work, w is in- 
dependent of L at high disorder and proportional to L 
at low disorder. If one assumes that w diverges at A c in 
the thermodynamic limit, w oc (A — A c )~ v , then one can 
determine A c and v from finite-size scaling collapses of 
the fingerwidth data. While this assumption was used in 
previous studies, our results with larger system sizes in- 
dicate that it is not justified. Finite-size scaling collapses 
become worse and worse as the range of L increases, and 
give estimates for A c that are clearly in the self-affinc 
regime. 
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FIG. 4. (a) Fingerwidths w calculated as a function of dis- 
order A for the simulation cell sizes indicated in (b). The 
dashed line gives an upper bound, obtained from the statisti- 
cal probability of singlets or pairs of unflippable s = — 1 spins. 
The solid line shows the external fingerwidth w e for L = 768. 
Panel (b) shows w e vs. A for all L. In both panels finger- 
widths were calculated in the plane normal to the nominal 
growth direction, z. 



The reason that w does not diverge at A c is quite sim- 
ple. Even in the self-affine regime there are small clusters 
of unflippable spins (Fig. ||) due to the tails in the dis- 
tribution of rji. These clusters are left behind by the 
advancing interface, and do not affect the morphology at 
long-length scales. However they do lead to a finite value 
of w. The limiting value of uu can be estimated from the 
probability for single isolated spins, pairs of spins, etc.. 
The upper bound for w from single spins and pairs is in- 
dicated by a dashed line in Fig. 0(a). This upper bound 
clearly inhibits divergence of the fingerwidth with L at 
any finite A. It is only about 300 at the value of A c 2.5 
determined below, and lower bounds would be obtained 
by considering larger clusters. 

We attempted to define improved fingerwidths by elim- 
inating all isolated spins, and then all pairs of isolated 
spins before determining the fingerwidth. However, this 
procedure was inefficient and did not converge rapidly. 
In the following we focus entirely on the morphology of 
the external interface, and thus eliminate surrounded un- 
dipped regions of all sizes. The scaling behavior in the 
self-affine regime is only associated with the external in- 
terface: Once the surrounded regions are left behind, 
they become irrelevant. In the self-similar regime we 
know that the external interface of a percolation clus- 
ter has the same fractal dimension as the entire cluster. 
Thus the external interface should give nearly the same 
fingerwidth as the cluster. 

The solid line in Figure [|(a) shows the fingerwidth cal- 
culated from the external interface, w e , as a function of A 
at L = 768. As expected, the unflipped regions have no 
effect on fingerwidth for A > 3, where the interface forms 
a fractal percolation pattern with narrow fingers. How- 
ever as A decreases towards A c « 2.5, the solid line rises 
sharply above the other data points. Figure [|(b) shows 
how the external fingerwidth changes with L. Finite-size 
scaling collapses of w e are discussed below. 



2. Interface roughness 

Self-affine interfaces are characteriz£d by the scaling 
properties of the interface roughness. Ela Due to our pe- 
riodic boundary conditions, the average direction of the 
external interface is normal to z. Its position is given by 
a height h(x, y) that may be multivalued. The roughness 
p(£) over a square region of side I in the x — y plane can 
be quantified by the root-mean-squared (rms) variation 
in h 



P (i) = ^((h-{h)r) e 



(3) 



where (h) is the average over a given square and {)e 
indicates an average oveE-||all square regions of side I. 
For a self-affine interface,^ p(£) ~ l a at large I. The 
roughness exponent a < 1 characterizes the degree of 
anisotropy, and would be unity for a self-similar fractal. 



4 



A = 2.1 




x 



FIG. 5. Cross-section through a system of size L = 96 for 
A = 2.1. The circles (open and solid) show all spins on the 
external interface. Solid circles indicate the single-valued in- 
terface obtained by taking the highest spin at each point in the 
(x,y) plane. Overhangs are evident in the regions where the 
solid and open circles separate. These form as the interface 
grows around unflippable clusters (crosses). As discussed in 
the text, unflippable clusters prevent the total fingerwidth w 
from diverging. However, the interface in this figure is clearly 
self-affine, and the largest unflippable cluster is much smaller 
than the system size. 

Figure [| shows a cross-section through a system of size 
L = 96 for A = 2.1, which is well into the self-affine 
regime. The external interface (circles) is a multi-valued 
function. In some regions there are overhangs where the 
interface extends over itself. These overhangs are neces- 
sary if the interface is to grow around small clusters of 
unflippable spins (crosses). The size of overhangs and of 
clusters of unflinped spins increases as A rises to A c . 

Previous workEH shows that overhangs can change the 
scaling of p(£) at small I. One way of highlighting their 
effect is to compare p(£) with the roughness pt(£) of the 
single- valued interface, h t (x,y), obtained by taking the 
top (highest) point on the external interface at each £ 
(closed circles in Fig. |]). Figure || shows log- log plots 
of both quantities vs. £ at the indicated values of A. 
For A < 2.2, results for the external (solid symbols) and 
single- valued (open symbols) interfaces converge at large 
I. The slope of the curves in the converged region is 
consistent with the roughness exponent-. a = 2/3 that 
is predicted from scaling^rgumentsEFa and observed 
in previous simulations. ErEu At small £ the dashed and 
full lines separate due to overhangs. The value of p t 
goes to zero at £ = 1, while the value of p goes to the 
rms variation in height above a single point in the (x, y) 
plane. The growing separation between dashed and full 
lines as A increases towards A c = 2.52 shows that the 
size of the overhangs increases. For A = 2.5 there is 



no convergence of the lines even at the largest £ and L 
we could study. The overhang size is one measure of a 
diverging length as A approaches A c from below, and is 
analyzed in following sections. 




FIG. 6. External interface roughness at the indicated val- 
ues of disorder A for L — 768. The full lines and symbols 
correspond to the total interface p, while the dashed lines 
and open symbols refer to the single- valued top interface pt. 
The straight solid line below the curves shows the slope cor- 
responding to a = | . The inset at the lower right shows that 
p and pt converge for A = 2.1 when a larger system size is 
used (L — 1260). There is some rounding of the curves at 
large £/L due to finite-size effects. 

3. Determining A c 

The data in Figures || and |^ give clear evidence of a 
morphological transition, and rough bounds on the value 
of A c . The saturation of the external fingerwidth with 
increasing L for A > 2.65 gives an upper bound for A c 
(Fig. ||(b)), while the merging of p and p t at large £ for 
A < 2.2 gives a lower bound (Fig. ^). In this subsec- 
tion we investigate other morphological attributes of the 
external interface that provide more accurate bounds for 
A c . In the next subsection we determine the exponent v 
through finite size scaling analysis of these quantities. 

A single-valued interface, hb(x,y), can also be defined 
by taking the bottom (lowest) value of the external in- 
terface h for each (x,y). Figure ^(a) shows the aver- 
age heights of the top and bottom interfaces, (h t ) and 
(/ib), as a function of A at various system sizes. In the 
self-similar regime, the fractal external interface extends 
throughout the entire height of the cell. The value of 
(ht) is a constant fraction of the system size and (hb) is 
of order of the fingerwidth. The two averages converge in 
the self-affine regime, where the difference between them, 
dh = (h t ) — (hb), is a measure of the height and abun- 
dance of overhangs. 
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FIG. 7. (a) Average heights of the single-valued interfaces 
(h t ) (dashed lines) and {hi) (solid lines) as a function of A 
for system sizes ranging from L = 48 (crosses) to L = 768 
solid triangles), following the symbol definitions in Figure 
(b) Average height difference dh = {ht — ht), normalized 
by system size, as a function of disorder. The crossing point 
where curves for all L intersect gives an estimate for A c . 



In Figure ^(b) we plot dh/L vs. A. As implied by the 
above discussion, this ratio vanishes at small A and rises 
to a constant fraction at large A. The increase becomes 
sharper with increasing L, and there is a clear crossing 
of all curves at A w 2.5. This means that for A below 
the crossing point dh/L decreases with L, while above 
the crossing point dh/L increases with system size. We 
conclude that the crossing point must coincide with A c . 

We have examined a variety of other quantities to con- 
firm that all give consistent values of A c and to minimize 
the error bars. Figure |8| shows results for two probabili- 
ties that are related to the global minimum, of each 
external interface. Data points connected by solid lines 
give the probability that h-/L is greater than 1/3. This 
probability is unity in the self-affinc limit and drops to 
zero in the self-similar regime where the fractal external 
interface extends all the way to the bottom of the system. 
The data points connected by dashed lines in Figure ^ 
give the probability that remains at the height of the 
initial seed plane. This probability is unity in the self- 
similar regime and drops to zero in the self-affine regime. 
Both probabilities exhibit sharper transitions from one 
to zero as L increases, and should become step-functions 
at A c in the limit L — > oc. Crossing points for the two 
probabilities in Fig. |[ the interface width in Fig. 0(b), 
and all other quantities that we examined are consistent 
with A c = 2.52 ±0.03. 

It is interesting to note that Fig. || shows a maxi- 



mum in H C (A) at A c . This is a reasonable result, given 
the difference in growth mechanisms for self-affine and 
self-similar regimes. In the self-affine regime, the inter- 
face must advance across the entire width of the system. 
Thus H c is sensitive to the regions that are hardest to 
flip, and rises with A. In the self-similar regime, the 
interface follows the path of least resistance. Since the 
number of spins that must be flipped (the percolation 
probability) is less than 1/2, increasing A makes it eas- 
ier to flip enough spins to span the system, and decreases 
H c . The rate of decrease in H c can be calculated exactly 
from the percolation probability in the large A limit, 
where spins are decorrelated.ErLl The asymptotic slope is: 
dH c /dA = -0.4907. 




FIG. 8. Probabilities related to the global minimum /i_ of 
the full external interface. The data points connected by solid 
lines give the probabilities that h- > L/3, while the data 
points connected by dashed lines give the probability that 
at least one point on the external interface remains at the 
initial seed plane. System sizes range from L — 48 (crosses) 
to L — 768 (solid triangles), following the symbol definitions 
in Figure ^| As L increases, the probabilities show a sharper 
change from zero to one. Quantities were averaged over an 
ensemble of at least 100 samples at each value of A and L. 



4- Finite-size scaling determination of v 



As in Sec. Ill A, we use finite-size scaling to determine 



the exponent v that describes the diverging correlation 
length £ at A c . The deviation from A c is measured by 
S = (A - A c )/A. We assume that £ - 8~ y , and that 
close to the critical disorder the only relevant lengths are 
£ and the system size L. Then dimensionless quantities 
like those shown in Figs. 171(b) and || can only depend 
on L/£, or equivalently LP U 5. When plotted against 
L 1 >"5, results for all system sizes should collapse onto a 
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universal scaling function. 

Figure § shows a scaling collapse for the data of Fig. 
0(b). The scaled interface widths dh/L for system sizes 
L = 48, 96, 192, 384 and 768 collapse well onto a univer- 
sal curve near L 1 ^ '5 — 0. As IL 1 / 1 ^! increases, the curves 
for small L begin to deviate from the others. These de- 
viations reflect corrections to scaling. They appear first 
at small L because these data points are for larger values 
of 8 than their counterparts at large L. 

Since the magnitude of corrections to scaling is not 
known, there is some uncertainty in determining the val- 
ues of A c and v. We found acceptable collapses for dh/L 
with A c = 2.52 ± 0.03 and u = 2.5 ± 0.3. Scaling col- 
lapses of other quantities related to the interface width, 
including p(l), the ratio of the lowest and highest points 
on the entire interface, and the probabilities shown in 
Figure ||, all gave consistent ranges of A c and v. We 
also considered the scaling variable S^^JA — A c )/A c , 



which gives different, and often larger ,113113 corrections to 
scaling. This led to a narrower region of scaling, but the 
same range of values for v and A c . 

The scaling behavior of the external fingerwidth is 
more complicated, because the distribution of finger- 
widths becomes bimodal in the self-affine regime. Most 
of the fingerwidths are essentially equal to the system size 
L. However, there is a significant fraction of very small 
fingerwidths from the region of width dh near the top of 
the interface (see Fig. ||). These make a disproportionate 
contribution to w e that does not obey the scaling ansatz. 
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FIG. 9. Finite size scaling collapse of the data in Figure [?] 
(b) for A c = 2.52 and v = 2.4. The symbols used for each L 
are the same as in Figure ^. 

The contribution of small fingerwidths decreases if one 
calculates higher moments of the fingerwidth. We define 



where n — 1 gives the mean width, n — 2 gives the rms 
width, etc.. We find a steady improvement in finite-size 
scaling collapses with increasing n. Figure [l0| shows that 
results for w e ^/L collapse onto a universal curve at large 
L with the same A c and v used in Fig. Best fits 
for v increased consistently from 2.0±0.2 at n = 1 to 
2.2 ±0.2 for n — 4, and the quality of the collapse showed 
progressive improvement. 

If w e is proportional to the diverging correlation 
length, then it must diverge with the same exponent. To 
check this, we examined the slope of plots of log 10 w en 
against log 10 S. The slopes were indeed consistent with 
values of v from finite-size scaling, although the uncer- 
tainties were somewhat larger. 



o 




en = Vy^T 



(4) 
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FIG. 10. Finite size scaling collapse of w e 4 using A c = 2.52 
and v = 2.4. The symbols used for each L are the same as in 
Figure ^ 



IV. SUMMARY AND DISCUSSIONS 

We have studied the zero-temperature phase diagram 
for interface growth in the 3D RFIM with a Gaussian 
distribution of random fields. There is a single transition 
from self-affine growth below A c to percolative growth 
at larger disorder. The tails of the Gaussian distribu- 
tion eliminate the faceted regime obseiafed in previous 3D 
studies of bounded field distributions B This means that 
the Gaussian RFIM provides a more realistic description 
of transitions in systems that do not have an underlying 
crystalline lattice, such as random porous media or amor- 
phous magnets. The phase diagram is also quite different 
from that for the 2D RFIM where Gaussian randomness 
suppresses the self-affine regime. 

The critical behavior at the onset of motion in the 
self-affine and self-similar regimes was analyzed using 
finite-size scaling. The critical exponent v>h that describes 
the diverging length scale as H — > H c was found to be 
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0.75 ± 0.02 for self-affine growth and 0.88 ± 0.02 for self- 
similar growth. These values are consistent with results 
for bounded distributions of disorderfl We also found 
the same roughness exponent a — 2/3 in the self-affine 
regime (see Fig. [|) . These results indicate that changing 
the form of the distribution of random fields does not 
change the universality class of the self-affine and self- 
similar growth regimes. 

The multi-critical point that separates self-affine and 
percolative growth was also analyzed. We found that the 
fingerwidth used in previous work does not diverge and 
can not be used to determine A c . Examination of the 
external interface revealed two lengths that did diverge 
at A c : The overhang size dh diverges as A increases to A c 
in the self-affine regime, and the external fingerwidth w e 
diverges as A decreases to A c in the percolative regime. 
Finite-size scaling collapses of these and other quantities 
gave consistent values for A c = 2.52 ± 0.03 and v = 
2.4±0.4. The error bars on these quantities are estimates 
of systematic uncertainties due to corrections to scaling. 

The value of v determined previouslyfj for a bounded 
distribution random fields, v = 3.0±0.5, is consistent 
with our result. However, this value was determined 
from the fingerwidth and is not reliable. Future work 
is needed to determine whether bounded and unbounded 
distributions are-,Hi| the same universality class. 

Perkovic et afitL3 have determined the critical behav- 
ior for the 3D Gaussian RFIM using a growth algorithm 
that appears to be in a different universality class. They 
analyzed the integrated avalanche size distribution oc- 
curring in one branch of a hysteresis loop (H increasing 
from — oo to oo). From the divergence as A decreased 
to A c they found numerical values for the critical disor- 
der, A c = 2.16 ± 0.03, and correlation length exponent, 
u= 1.43 ±0.18. 

The discrepancy between our results and those of 
Perkovic et al. seems to result from a crucial differ- 
ence in our growth algorithms. They allowed any spin- 
flip that lowered the energy, while we only allowed spins 
on the interface to flip. The exchange coupling between 
neighbors dominates in the low disorder limit, and spins 
are very unlikehtto flip unless they are on the interface. 
Nowak et al.li3tL3 found that the difference between the 
two algorithms was negligible in simulations of low dis- 
order growth with uniform disorder. The two algorithms 
should also yield the same percolating cluster in the high 
disorder limit, where interactions between neighbors be- 
come irrelevant. The problem maps onto ordinary perco- 
lation, and the order in which spins are flipped becomes 
irrelevant. Near A c these arguments break down, and the 
algorithms may give different results. At intermediate 
disorder, the exchange coupling is weak enough to allow 
clusters to flip ahead of the interface, and correlations are 
important enough that these flipped clusters can aid the 
advance of the approaching interface. One expects that 
both H c and A c will be lowered by the advance clusters 
in Perkovic et aZ.'s model, and this is consistent with the 
numerical results. 
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